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Abstract 

Spatial resolution of most imaging devices is fundamentally restricted by diffraction. This limita- 
tion is manifested in the loss of high spatial frequency information contained in evanescent waves. 
As a result, conventional far-field optics yields no information about an object's subwavelength 
features. Here we propose a novel approach to recovering evanescent waves in the far field, thereby 
enabling subwavelength-resolved imaging and spatial spectroscopy. Our approach relies on shifting 
the frequency and the wave vector of near-field components via scattering on acoustic plionons. 
This process effectively removes the spatial frequency cut-off for unambiguous far field detection. A 
straightforward extension of this technique, which we call spatial Fourier transform infrared spec- 
troscopy, allows to preserve phase information, making it possible to perform 3D subwavelength 
imaging. We discuss the implementation of such a system in the mid-IR and THz bands, with 
possible extension to other spectral regions. 



Optical microscopy is widely used for probing the nanoscale, with particularly important 
applications in the life sciences. It is fast, noninvasive, and can be coupled with fluorescence 
and spectroscopy studies to uncover the structure, composition, and dynamics of biological 
samples. However, due to the diffraction limit much of the important structural informa- 
tion about the object under study becomes lost [1]. We illustrate the extent of this loss by 
comparing the image [Figjlj^b)] of a subwavelength optical target [FigQa)] with one that 
would be produced by an ideal diffraction-limited optical system (which could be imple- 
mented using e.g. synthetic aperture techniques [21 E])- It is apparent the details smaller 
than ~ A/2 cannot be resolved. Because of this wavelength dependence, the resolution 
of the image particularly suffers in the mid-IR and THz spectral bands. This creates a 
significant stumbling block for studying microscopic chemical morphology. For instance, 
mid-infrared spectroscopy - a crucial tool for chemical identification - cannot be used for 
structure-specific analysis of biological samples due to the inability to resolve small cells, 
bacteria, and cellular subunits. 

There exist certain techniques that can improve spatial resolution of optical imaging be- 
yond the diffraction limit. One category of methods involves scanning probes |1H6], which 
rely on the detection of evanescent waves in the near field of the sample. Nonlinear processes 
form another set of approaches; these include coherent anti-Stokes Raman spectroscopy 
(CARS) [TIIH], as well as RESOLFT techniques (which utilize saturable fluorescence transi- 
tions) [0]. All the above solutions typically involve complicated and expensive experimental 
setups. This limits their widespread availability, as well as their potential for creating com- 
pact integrated detection and/or imaging devices. In addition, some setups may be difficult 
to adapt for operation in mid-IR or THz. 

In the past few years, there emerged several intriguing approaches towards subwavelength 
imaging that involve using metamaterial-based devices, which promise to be low-cost and 
readily adaptable for a variety of spectral regions. Such devices work by recovering the 
evanescent waves and the concomitant information about the object's fine structure. A 
negative index "superlens" [lOJ, for instance, can drastically amplify the evanescent com- 
ponents of the spatial spectrum while transmitting the propagating waves. However, the 
superlens exhibits exponential sensitivity to losses, limiting practical devices to near-field 
operation [TT]. A more recent set of experiments involved subwavelength scatterers placed 
in the near field of a source. By diffracting off the scatterers, evanescent waves convert into 
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propagating waves which could then be used to gather information about the near-field spec- 
trum [12] or to attain subwavelength focusing [13]. However, in this far- field superlens [12], 
the diffracted evanescent waves mix with the existing propagating spectrum leading to the 
loss of information about the target. 

In the present Letter, we propose an alternative approach to far-field spatially resolved 
spectroscopy, based on a device that converts mid-IR or THz evanescent waves to propagat- 
ing waves via scattering on acoustic phonons. These scattered and frequency-shifted waves 
can be easily decoupled from the existing propagating spectrum that forms the regular 
diffraction-limited image, and with minimal processing can be used to distinguish subwave- 
length features. Furthermore, the ability to dynamically tune the period of the acoustic 
grating makes this system robust, flexible, and capable of operating over a wide range of 
spatial frequencies. 

The proposed super-resolution Fourier microscope/sensing system is shown in Fig. |2} 
The object is placed in the near fleld of an acousto-optic modulator (AOM) and illuminated 
with a plane wave from a mid-IR or THz source. Waves scattered from the object strike 
the phonon grating set up in the AOM by a running acoustic wave with frequency Q. Due 
to scattering on the phonons, the transverse wave vector kx of the incident radiation is 
shifted by integer multiples the phonon wave vector q, while its corresponding frequency is 
shifted by integer multiples of Q. For a sufficiently large q, the evanescent components of 
the object's spatial spectrum (|A;^| > u/c) can be scattered into the propagating waves with 
\k'^\ = \kx — q\ < oj/c. The various spatial frequency components can be measured using a 
Fourier optics setup (e.g. a lens with a detector array in its focal plane). 

The amplitudes of the frequency-shifted waves can be obtained using the standard meth- 
ods of acuostooptics [T3] (cf. Methods). We define A± = f^Ei^^kxTQ), = ^o-Sin(fcx) (with 
linear coefficients to, t^ describing the generation of phonon-scattered and/or device trans- 
mission characteristics) and assume Ai ^ Aq ^ A"^, where Ai is the detected amplitude of 
the illuminating wave Aie^^°^ . Averaging out the signal over the finite detector aperture and 
subtracting the background (which can be done electronically), we can write the intensity 
detected by the system of Fig. |2] as 

/out = |ioP + 2 + |i,i+|2y^'cos(f]t + 7). (1) 

The two terms in this equation can be decoupled using standard techniques: the DC term 
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is isolated with the aid of a low-pass filter, while the term oscillating at the acoustic fre- 
quency Q is recoverable using lock- in detection. For any given k^, this second term contains 
contributions from both = f^E^T^ik^ — q) and y4„ = t~ E-^T^ik^ + q). Although the cou- 
pling between these two quantities, together with the lack of phase information, makes it 
difficult to recover the spatial spectrum, the information collected can be used in detecting 
subwavelength morphological changes between different samples. 

To illustrate this, we utilize Eq. ([T| to perform a comparison between the standard 
optical target of Fig. [T]and a modified target, where the label of every 6th line group has 
been randomly replaced. The first replacement corresponds to the last resolvable line group 
(A/2.5 line separation); the subsequent replacements correspond to halving the size of the 
line groups (A/5,. . ., A/40). We assume the measurement is performed by selecting an 
element of a photodetector array in the observation plane and using two orthogonal acoustic 
transducers to scan the acoustic wavevector within the range q^^y G [— 25(X'/c, 25(X'/c]. 

It should be emphasized that any method that relies on digital processing of raw data 
can suffer form rapid - sometimes exponential [TT] - accumulation of noise. To address this 
potential issue, in our computations we add a normally-distributed random term to the AC 
amplitude of Eq. ([T]) in order to simulate noise in the system. Because SNR is expected 
to be lowest for maximum values of the acoustic wavevector g, we consider SNR=10 for 
q = 25u/c |23]. Assuming a practical 20x20 element photodetector array, we compute the 
signal given by Eq. ([T]) for the standard target, as well as the modified target [Fig. |3]^a)]. 
Fig. |3|^b) shows the result of subtracting the two datasets and performing an inverse Fourier 
transform, with the resulting plot superimposed onto the modified optical target. Evidently, 
every change in the original image is manifested in this difference diagram. Furthermore, it 
is largely localized in the vicinity of the actual changed pixels. It is possible to discern the 
difference signal even from the A/40 line group label. 

The ability to distinguish between fine spatial features of optical targets makes the system 
described above uniquely suited for identifying objects based on their subwavelength spatial 
features. As a result, it may find many applications in fingerprinting and/or detection of 
chemical and biological structures. 

Furthermore, a straightforward modification of this setup not only allows to measure 
the "downshifted" A_ component directly, but also provides a method for retrieving phase 
information, making it possible to perform phase-contrast microscopy, as well as 3D imaging 
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on subwavelength scales. 

To this end, a portion of the illuminating radiation is shifted in frequency by Qi, using a 
second AOM. Unlike the modulator that interacts with light scattered from the sample in 
the Raman-Nath regime [15] , this second AOM utilizes an appropriately oriented and longer 
cell to produce Bragg scattering. This results in a strong optical signal at frequency u + Qb, 
\Ai,\ exp[i{kb-r — {u + Qb)t)], which is projected onto the detector [see Fig.|4|a)]. Interference 
between the two optical signals produces beat note photocurrents with frequencies Q, Qb, 
ilb + ^b- ^■ 



loutiK) = Ei exp{ikoz) + Ab exp{ik ■ r) 

+ [A~ exp(zfit) + exp{-int) + Aq] exp(iA; ■ r) 

= . . . + 2|i~ib| cos[{nb + n)t + A$^] 

+ 2\A+Ab\ cos[(^]b - n)t + A$+] + . . . , (2) 

where = {kb — k) ■ r — (p^ is the phase difference between the signal from the Bragg cell, 
|y4b| exp(i/cf,-r), and the Raman-Nath-scattered signal exp(z/c -r) = exp[i(0^ + /c-r)]. 

Of special interest is the component at frequency Vl + Vlb-, which carries the high spatial 
frequency information contained in its modulus and its phase A$~ ~ {k^j. — kx)x—(f)~. Both of 
these quantities can be retrieved using lock-in techniques. To produce the lock-in reference, 
the RF signals driving the two acoustic cells can be mixed using a nonlinear element (e.g. 
a diode) and appropriately filtered to produce the sum frequency. As a result, complete 
information can be obtained about the complex high spatial frequency Fourier component 
A", from which it is straightforward to deduce the field E\-a{kx + q). By collecting data 
from multiple CCD pixels, as well as by varying the acoustic wave vector g, information 
can be collected about the spatial spectrum of the object. The data can then be digitally 
processed to produce a spatial-domain image containing subwavelength details, as well as 
phase contrast. 

Because the Bragg-shifted signal we use to decouple the Aj^ and A_ terms serves as a 
reference needed to record phase information, and because the image is reconstructed digi- 
tally, our technique bears some similarities with digital Fourier holography (DFH) [T614T8] . 
However, our method contains several key enhancements over DFH. In conventional holog- 
raphy, care has to be taken to isolate the target signal both in real and Fourier space. This 



constrains reference wave geometry, translating into limitation on the field of view, as well 
as maximum attainable resolution. The requirement that the CCD pixel spacing must allow 
for imaging the reference wave fringes further limits the resolution. By virtue of frequency- 
shifting the signal, it is possible to isolate the interference term of interest. Furthermore, 
since the spatial spectrum measurements are performed not only by selecting different CCD 
pixels, but also by scanning the acoustic wavevector, the limitations of CCD's physical 
spatial frequency bandwidth (introduced by pixel granularity) [17] can be circumvented. 

We simulate the performance of the system by first using Eq. ^ to compute the response 
of the system to a calibration signal having unit amplitude for all spatial frequencies. In 
practice, such calibration signal might be generated by placing a point source in the vicinity 
of the AOM. Eq. ^ also provides the effective amplitude and phase transfer functions that 
allow to determine the detected signal for a given input field distribution. Gaussian noise 
is added to simulate spurious signals in the system. The input signal can then be obtained 
by dividing out the calibration quantities. In Fig. |4](b) we plot the simulated retrieved field 
magnitude. Zooming in on the central part of the test patten (figure inset) it is evident 
that every line group is distinctly resolved. Furthermore, because the phase information is 
preserved, the full 3D information about the target is collected. 

There exist many possible ways to enhance the functionality and the performance of the 
proposed devices. For instance, sensitivity of the device may be improved with a subwave- 
length layer of highly doped semiconductor at the front AOM facet. When the dielectric 
constant of this layer is equal to —1, the evanescent fields are strongly enhanced due to 
resonant coupling to surface plasmons in the doped layer (a phenomenon known as "poor- 
man's superlensing" [ID]), leading to better SNR at the detector. Another possible way to 
improve the scattering efficiency of evanescent waves is placing the sample directly in the 
path of an acoustic wave, for instance, by running the wave through a microchannel con- 
taining objects to be studied. This approach may find many applications in novel integrated 
biological/chemical detection devices. Finally, we comment on the possibility of extend- 
ing the proposed approach to frequencies other than the mid-IR and THz bands discussed 
here. While implementing the system for lower frequencies is essentially trivial, near-IR and 
optical frequencies pose a challenge. Acoustic phonon energies in practical devices do not 
approach the values necessary to produce a substantial wave vector shift in these spectral 
bands. However, the required shift in spatial and temporal frequencies can in principle be 
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attained by replacing the acoustooptic medium with a nanostructured periodically moving 
grating. 

In conclusion, we have proposed a system that enables detection of sub-diffraction-limited 
spatial spectrum components in the far field by utilizing scattering from an acoustic grat- 
ing. This process works whenever the spatial frequencies of the object are comparable in 
scale to the acoustic wave vector. In its simplest implementation, the system could aid in 
"fingerprinting" of samples based on their subwavelength spatial features. With the use of 
an additional Bragg-shifted reference signal, it is also possible to recover the phase of the 
original optical signal. This technique, the spatial Fourier transform infrared spectroscopy, 
allows to perform subwavelength-resolved 3D imaging. In addition, the proposed approach 
has the potential to greatly enhance the specificity of mid-IR and THz spectroscopy. 

Details of the calculations 

We start by considering a rectangular sound column (i.e. planar acoustic wavefronts 
propagating in the x direction) interacting with a spectrum of incident plane waves in a 
dielectric medium. We neglect the diffraction of the sound field and assume weak interaction. 
Due to photoelastic effect, the sound field produces a sinusoidal modulation of the dielectric 
permittivity [15j, 



which corresponds to a spatiotemporal volume grating. We may write the general form of 
the field inside the grating as a sum over the discrete diffracted orders. 



j 

The scattered plane wave components Aj{z) are then governed by the Raman-Nath equa- 
tions [n]. 




(3) 




(4) 



A';{z) + klA,{z) = 



(5) 



where k^. 




1 1/2 

{kx+jqy , and 
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The amplitude of the j*^ diffracted order, Aj, is proportional to (Ae)-^, with Ae/e ^ 1, 
allowing to ignore higher order terms (j > 2). We can, furthermore, conclude that the 
amount of energy scattered into the shifted waves is small, thereby permitting to neglect 
the variation of O*'^ diffracted order Aq (the undepleted pump approximation) [15]. We 
note that for propagating waves, this conclusion is valid insofar as there exists no Bragg 
matching between the incident and diffracted waves. Since the phonon wave vector g is a 
tunable parameter in our model, it is always possible to pick a range of q values to ensure 
minimal energy loss in the incident wave. For the evanescent waves, the undepleted pump 
approximation is justified by the small interaction length. 

Keeping terms up to first order in Eq. Q, we obtain: 

A'^{z) + klA^{z) = -^^uj]Ao{z). (6) 

The scattering amplitudes of "upshifted" and "downshifted" waves can be obtained from 
this expression. Since the input field at spatial frequencies (kxTl) contributes to the output 
field at spatial frequency k^:, we may write: 

Eout{kx) = A_ exp{iQt) + A^exp{—iflt) + Aq exp[i{kxX — ut)], (7) 

with A± = t^Ei,^{kx =F ^0 = ^o^o- Scattering coefficients are given by 



2 \c) q{qT'2k^) "^k. 



where kf = a/ e{uj/cY — [k^ =F qYi kz = a/ e{uj/cY — k^, and to is a transmission coefficient. 
We note that for g ^ 0, as well as g ~ i-2kx, the perturbative treatment of Eq. ^ breaks 
down due to the onset of Bragg-matching. 

From Eq. ^ we can estimate the diffraction efficiency of high spatial frequency input 
signal Ei^{k]^) as 

uVc Ae 

' ' 2kf n(l + n) ' ^ ^ 

with n = a/c (the refractive index of the acoustic medium), and Ae oc VF, the flux of 
acoustic energy per unit area. 

In our computations, we assume the operating wavelength of 10 fim with germanium as 
the acoustic medium. We take Ae = 10^^ and restrict the magnitude of the acoustic wave 
vector q to 25u/c. To obtain Ae = 10^'^ the ultrasonic fiuence of 33 W/cm^ is required. 
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Since for high spatial frequencies fc™ ~ q, acoustic driving frequencies up to 8.75 GHz are 
required to retrieve fc™ ~ 25uj/c. These parameters are within reach of modern uhrasonic 
transducers p^, as well as surface acoustic wave devices [20]. the system depends on the 
characteristics of the The signals are generated by sampling 10® points in spatial frequency 
space. Assuming 400 points can be sampled at once with a detector array, it is necessary 
to use a focusing lens with N.A. ^ 0.6 and adjust the acoustic wave vector components 
(lx,y £ [— 25,25] cj/c sequentially, effectively scanning the low N.A. system over a larger 
spatial frequency spectrum [31 EI]- At KHz readout rates, short (< 1 s) acquisition times 
can be obtained with this setup. 

We simulate the measurements by taking the signal indicated by Eqs. ([T]) and ^ and 
introducing additive Gaussian noise. Assuming a shot-noise limited long-wavelength IR 
detector with a typical detectivity D* ~ lO^cmv^Hz/W [22] we find that ^ 20 — 30 mW of 
illuminating optical power is needed to obtain SNR of 250//cJP used in our calculations. 

Finally, we note that a potentially substantial source of noise in the proposed system 
is the detection of zero-order (undiffracted) illumination at the shifted frequencies due to 
the finite source linewidth. In the measurement region, the Lorentzian lineshape of the 
source takes the same functional dependence as the diffracted signal of Eq. (|9]). Thus, 
the zero-order illumination simply adds a constant measurement background that may be 
subtracted. The noise floor is effectively raised by a factor '^/^ n{i+n) ' where = v/\ is 
the acoustic frequency that yields q = u/c, and AQ is the source linewidth. Assuming the 
source is a frequency-stabilized quantum cascade laser with a ~ 15 KHz linewidth [23], this 
increase factor is not significant (< 2). 

The work was partially supported by ARO MURI and NSF MIRTHE. 
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FIG. 1: Computed effects of the diffraction limit on an optical test target (assuming coherent 
illumination). The target consists of labeled sets of line groups, which halve in size every six 
elements. The last resolvable set of lines (group 5, left column) corresponds to 1.5 lines/A. 
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FIG. 2: Schematics of the proposed system. 
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FIG. 3: (a) Optical test target and its modified version (inset). In the modified target, the "5" 
label of every column has been replaced by another digit, (b) Computed output of the system 
in the presence of noise (shown in grayscale) assuming a realistic, noisy detector with 400 active 
photocells. The modified optical target is superimposed for illustration purposes. The output of 
the system clearly identifies the location of every modified digit, even for regions far below the 
diffraction limit. 
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FIG. 4: (a) Schematics of the proposed system. Note the Bragg-shifted reference beam that aids 
in providing phase information, (b) The computed output of the system with optical test target 
as the object in the presence of noise. 
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